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I compute the leading contribution to the ground state Renyi entropy Sa for a region of linear 
size L in a Fermi liquid. The result contains a universal boundary law violating term simply related 
the more familiar entanglement entropy. I also obtain a universal crossover function that smoothly 
interpolates between the zero temperature result and the ordinary thermal Renyi entropy of a Fermi 
liquid. Formulas for the entanglement entropy of more complicated regions, including non-convex 
and disconnected regions, are obtained from the conformal field theory formulation of Fermi surface 
dynamics. These results permit an evaluation of the quantum mutual information between different 
regions in a Fermi liquid. I also study the number fluctuations in a Fermi liquid. Taken together, 
these results give a reasonably complete characterization of the low energy quantum information 
content of Fermi liquids. 
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I. INTRODUCTION 



Recently, motivated by studies of entanglement in 
quantum systems, there has been a profitable exchange of 
ideas between the quantum information and condensed 
matter communities. Much of this exchange has focused 
on the physics of entanglement entropy, defined as the 
von Neumann entropy of the reduced density matrix of 
a subsystem. Entanglement entropy is the basic quan- 
tity that provides a characterization of the entanglement 
properties of many body quantum statesii^. Insights 
gained from the study of entanglement led to a new class 
of variational many body states collectively known as 
tensor network states^"—. In addition, entanglement en- 
tropy can characterize interesting topological and critical 
phases of matter. Conformal field theories in 1-t-l dimen- 
sions are partially characterized by the leading universal 
term in the their entanglement entropy^^. On the other 
hand, most phases in more than one spatial dimension 
have an entanglement entropy that scales as the bound- 
ary of the spatial region considered^il. This boundary 
term is non-universal and typically forces us to study 
sub-leading terms in the entanglement entropy to find 
universal quantities characteristic of a particular phase. 

There exists a class of systems in higher dimensions 
that violate this boundary law for the entanglement 
entropy. This class includes free fermions^"— , Fermi 
liquidai^, Weyl fermions in a magnetic field'^^, and more 
exotic systems including critical Fermi surfaces and Bose 
metals. The apparent unifying theme is the existence 
of many gapless one dimensional degrees of freedom in 
all these systems. Fermi liquids with a codimension 
one Fermi surface have long been known to be equiva- 
lent to a set of nearly decoupled one dimensional gapless 
modesii"— . The system of free Weyl fermions in a back- 
ground magnetic field is also explicitly 1-1-1 dimensional 
due to the presence of zero modes propagating along the 
magnetic field lines. Even critical Fermi surfaces behave 
thermodynamically like a collection of 1 -I- 1 dimensional 
critical fields with critical exponent z/ 7^ 1?°. In all these 
cases we may think of one dimensional theories patched 
together to form a higher dimensional theory, but this 



picture is most concrete for the case of Fermi liquids. 

In this paper I use the one dimensional framework to 
compute the Renyi entropy, defined below, for interacting 
fermions in a Fermi liquid state. This quantity is inter- 
esting because it gives in principle complete knowledge 
of the spectrum of the reduced density matrix, although 
I will only compute it in the low energy limit. I also de- 
scribe the form of the entanglement entropy of more gen- 
eral regions (non-convex or even disjoint) using the same 
one dimensional formulation. This permits a calculation 
of the mutual information between two distant regions in 
a Fermi liquid. All of these information theoretic quan- 
tities turn out to provide direct access to the geometry 
of the interacting Fermi surface. They are universal in 
the sense that they depend only on the geometry of the 
interacting Fermi surface and not on any other details 
of the Fermi liquid fixed point. Additionally, I compute 
the fermion number fluctuations in a region of size L in a 
Fermi liquid and give general finite temperature crossover 
forms for both the entropy and the number fluctuations. 
Together these results provide a very complete charac- 
terization of the quantum information theoretic content 
of Fermi liquids in terms of the geometry of low energy 
excitations. 



II. ONE DIMENSIONAL FRAMEWORK 

The low energy physics of a system of free fermions is 
given in terms of the geometry of the free Fermi surface. 
For example, the heat capacity of such a free fermion sys- 
tem is simply proportional to the density of states times 
temperature, and the density of states may be written as 
an integral of l/vp over the Fermi surface. Each patch of 
the Fermi surface contributes to physical quantities like 
a chiral 1-1-1 dimensional free fermion. 

Including interactions is possible at the free Fermi fixed 
point because phase space restrictions reduce the effects 
of interactions to certain forward scattering terms labeled 
by an infinite set of Landau parameters. Each Landau 
parameter corresponds to a single exactly marginal de- 
formation of the free Fermi fixed point which preserves 



all scaling dimensions. These forward scattering interac- 
tions preserve the fermion number on each patch of the 
Fermi surface, and the Fermi liquid has a very large sym- 
metry group: t/(l)°°— 12^. This one dimensional frame- 
work permits computation of the usual physical observ- 
ables of Fermi liquids, observables like heat capacity and 
compressibility. In addition, it provides simple access 
to many of the anomalous entanglement and fluctuation 
properties of Fermi liquids. 

These anomalous properties characterize the reduced 
density matrix of a region of linear size L inside the Fermi 
liquid, and they represent non-traditional global observ- 
ables as opposed to local correlations functions. These 
observables include entanglement entropy, defined as the 
von Neumann entropy of the reduced density matrix of 
the region, and fermion number fluctuations, defined as 
the variance of the fermion density integrated over the 
region. The entanglement entropy can be generalized to 
the Renyi entropy defined as 
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The von Neumann entropy S'vn — ^^"^{pLlnpL) is re- 
covered from Sa in the limit a — ;■ 1. Fermionic sys- 
tems without a codimension one Fermi surface and most 
bosonic systems have a Renyi entropy which scales with 
the boundary L'^~^ of the region considered. This leading 
scaling behavior is sensitive to the cutoff of the low en- 
ergy effective theory and does not define a universal low 
energy observable. Other sub-leading terms in the entan- 
glement entropy may provide universal numbers charac- 
terizing different phases, but these terms are in general 
hard to calculate and interpret. 

Fluctuations of conserved quantities also behave sim- 
ilarly to the entanglement entropy. Consider a system 
of fermions with a conserved fermion number obtained 



by integrating the density n{x) over the entire system. 
Given access to only a subregion of linear size L we may 
ask about the observable Nl = JxeL-^ n{x). The average 
(iVi) — tr(pi,iVi) is simply the average density times 
the volume of the subregion for a translation invariant 
system. The fluctuations AA^I = {{Nl - {NL)f) = 
tr(pL(^L — (^l))^) are generically non-zero and typically 
scale as the boundary of the subregion. However, there 
are exceptions to this scaling, for example, it receives a 
logarithmic correction, like the entanglement entropy, for 
critical one dimensional systems. Also, symmetry break- 
ing states have fluctuations in conserved quantities that 
scale as the volume of the subregion; these fluctuations 
are due to the zero mode of the order parameter. 

Note that the presence of fluctuations in fermion num- 
ber within a subregion is not in conflict with a fixed to- 
tal fermion number for the entire system as we always 
study small subsystems of a larger system. Like the en- 
tanglement entropy of Fermi liquids, the number fluc- 
tuations in a metal with codimension one Fermi surface 
are anomalously large. One finds that /S.N\ ~ L'^^^ Ini 
with a multiplicative logarithmic correctio n^'^^i^^'^"^ . We 
will later determine the precise form of this leading log- 
arithmic term. 

These results are not restricted to zero temperature; 
indeed, there are universal crossover functions that de- 
termine how the entanglement entropy and ground state 
number fluctuations go over into the corresponding ther- 
mal quantities as the temperature is raised. For example, 
consider a circular region in a two dimensional spherically 
symmetric Fermi liquid. The reduced density matrix for 
this region displays crossover behavior at finite tempera- 
ture from the zero temperature anomalous entanglement 
regime to the finite temperature thermal regime. This 
crossover is captured by a universal scaling function 
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where Lcs is the chordal distance across the circle. This 
crossover function depends only on the geometry of the 
Fermi surface and the real space region; the physical 
Fermi velocity vp appears simply as a unit conversion. 
It is the only trace of interactions in this remarkably uni- 
versal formula. We will discuss and generalize this for- 
mula later in the paper. A similar formula exists for the 
number fluctuations which we will derive below. 



III. RENYI ENTROPY 

I turn flrst to the calculation of the Renyi entropy of 
the Fermi liquid. Like the entanglement entropy, the 
Renyi entropy is a cutoff sensitive quantity. It typically 
satisfles a boundary law, but we will flnd that Fermi liq- 
uids have a universal boundary law violating piece tied 
to the geometry of the Fermi surface. The Renyi entropy 
is of interest for several reasons. It is often easier to com- 
pute, both in field theory and numerically, for example, 
by quantum Monte Carlo2^. It can also be used to charac- 
terize phases of matter, for example in topological phases, 
conformal field theories, and holographic theories^ i^^'^^ . 



From a quantum information perspective, the notion of 
single copy entanglement^, given by the limit a ^' oo of 
Sa, determines the maximal deterministically distillable 
entanglement from the bipartition of a single copy of a 
quantum state. The entanglement entropy only gives the 
distillable entanglement in the unphysical limit of many 
copies of the quantum state. I find that the single copy 
entanglement is, for large L, simply one half the entangle- 
ment entropy, similar to the one dimensional conformal 
result^^. We will now use the one dimensional framework 
to compute the Renyi entropy of a Fermi liquid. 

Restricting first to the non-interacting case, each patch 
on the Fermi surface is equivalent to the chiral half of one 
dimensional free fermion. The Renyi entropy Sa of such a 
one dimensional chiral fermion is known from conformal 
field theory. In fact, it is simply proportional to the en- 
tanglement entropy with an a dependent prefactor. The 
precise result is 
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where L is the length of the one dimensional interval 
and e is a short distance cutoff. Returning to the Fermi 
surface problem, each patch contributes such a term to 



the entanglement entropy with cl — 1 and c/j = 0. The 
resulting Renyi entropy of the Fermi liquid is remarkably 
simple; the universal part controlled by the low energy 
theory is simply 
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where again. Si is the usual entanglement entropy in two 
dimensions given by 



Si = S = 
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This isn't the end of the story. Using the finite tem- 
perature crossover form for the von Neumann entropy of 
the reduced density matrix we obtain a finite tempera- 
ture crossover form for the Renyi entropy. The result for 
two dimensions is 
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and in any spatial dimension with a codimension one 
Fermi surface 
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I want to point out that since the Fermi surface system 
has a cutoff, this formula cannot truly be valid for all 
a. For example, the a ^ limit of the Renyi entropy 
always produces the Schmidt rank of the reduced density 
matrix under study. However, the a — ^ limit of the 
conformal field theory result diverges corresponding to 
the statement that a theory which is truly conformal to 
arbitrarily high energy must have an infinite number of 
local degrees of freedom. At finite temperature we must 
keep T/a < Tp or a > T /Tp in order not to probe 
high energy physics, and similarly, at zero temperature 
we must keep a > l/{kpL) to avoid the infiucncc of non- 
universal physics away from the Fermi surface. 



acting Fermi surface in this less symmetric systems con- 
tinues to control the Renyi entropy. Second, the Fermi 
velocity must be replaced by the renormalized physical 
Fermi velocity. 

There are several checks of this formula. First, it cor- 
rectly reproduces the one dimensional result that must 
occur for non-interacting systems with nested Fermi sur- 
faces. Second, it reproduces the finite temperature Renyi 
entropy of the free Fermi gas; this follows by a direct com- 
putation. Third, it reproduces the thermal Renyi entropy 
of an interacting Fermi liquid. At finite temperature the 
Renyi entropy is 



Now at last restoring the interactions, the complete fi- 
nite temperature crossover function remains correct even 
in the presence of interactions because of the nature of 
the Fermi liquid fixed point (more properly, fixed man- 
ifold). The counting of low energy degrees of freedom 
remains the same. There are two possible modifications 
of this formula when including interactions. First, the 
interactions may change the geometry of the low energy 
Fermi surface. Spherical symmetry can prevent such a 
modification, but in a solid state metal the crystal lattice 
breaks the rotational symmetry. Nevertheless, the inter- 
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where Z(/?) is the partition function. The partition func- 
tion of the Fermi liquid is —InZ^T for small T com- 
pared to Ep^ and this permits us to write Z{a/3) — 
ZiJiY ■ The final result is that the Renyi entropy for an 
interacting Fermi liquid is the same universal function of 
a times the thermal entropy as predicted by the formula 
above. 



IV. ENTANGLEMENT ENTROPY OF 
DISJOINT REGIONS 

Having described in detail the leading behavior of the 
Renyi entropy for a single convex region, I now turn 
to the problem of describing the entanglement entropy 
for more complicated subregion geometries. Perhaps the 
most important motivation for this study is the calcula- 
tion of the quantum mutual information I{A, B) between 
two regions A and B. The mutual information is 



I{A,B)^Sa + Sb-Saub, 



(9) 



and this definition requires knowledge of the entangle- 
ment entropy of disjoint regions to compute. We will see 
that the entropy of disjoint regions can be computed in a 
manner similar to that of convex regions via a more gen- 
eral formula that can handle arbitrary region geometry. 
The basic strategy remains the same: express all quanti- 
ties in terms of sums over one dimensional modes on the 
Fermi surface. Because these one dimensional modes are 
conformal and, in fact, essentially free fermions, much 
can be calculated. In particular, the main information 
we need is the entanglement entropy of a one dimen- 



sional conformal field theory when the one dimensional 
subregion consists of multiple regions. 

This quantity, the entanglement entropy of multiple 
regions in a conformal field theory, is not a simple quan- 
tity to compute in general. Despite some early claims 
in the literature it is not as universal as the single re- 
gion entanglement entropy which depends only on the 
central charge of the conformal field theory. The multi- 
interval entanglement entropy depends on the entire op- 
erator content of the theory as it requires evaluating 
higher point correlation functions of twist fields. How- 
ever, the relative simplicity of the Fermi liquid is a boon: 
the multiple interval result for a one dimensional free 
fermion is known. We will see shortly how this informa- 
tion permits calculation of the entanglement entropy of 
arbitrary regions in a Fermi liquid. 

Let us briefiy recall the one dimensional result. We 
consider a subsystem consisting of several disjoint inter- 
vals labeled [ai,bi\ i = l,...,m with m the number of 
intervals. As we said before, the entanglement entropy 
of such a subregion is in general complicated, but for free 
fermions the answer is remarkably simple. The result is 
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i<3 



with the implicit presence of cutoffs in the logarithms is 
understood. This formula reduces to the usual result in 
1 dimension in the case of a single interval. Note that 
when computing the mutual information, the number of 
logs being added and subtracted is equal giving a result 
independent of the cutoff. 

How do we use this formula to compute the entangle- 
ment entropy of a Fermi liquid in an arbitrary geome- 
try? The procedure is illustrated in Fig. 1 and goes as 
follows. First, we choose a point on the Fermi surface; 
this point defines a direction in space via the local Fermi 
velocity. Second, we choose a point on the real space 
boundary and draw a line through this point with di- 
rection determined by the Fermi velocity from step one. 
Third, we determine the intersections of this line with 
the boundary of the real space region; these intersec- 
tion points are grouped into effectively one dimensional 
intervals. The contribution to the entropy of such a con- 
figuration is given by a slightly modified version of the 
one dimensional result above pUj) . To each logarithm in 
the one dimensional result we append geometrical "flux 
factors" \nx ■ nk\ times a differential area element dAx 
where n^ is the direction determined by the Fermi veloc- 
ity and Ux is a real space normal vector. In detail, each 



logarithmic term of the form In |a: — j/| is multiplied by a 
factor of ^"°''"'° — ^+l"»-"fc — y^ where n^ and Uy are the 
local real space boundary normals. As illustrated in Fig. 
2, the flux factors and differential area elements satisfy 
\nx ■ nk\dAx — \ny ■ nk\dAy. Finally, the contribution 
of this configuration of intervals must be divided by the 
number intersection points to avoid over-counting and 
integrated over all such points on the real space bound- 
ary and the Fermi surface. This geometrical calculation 
gives the leading boundary law violating contribution to 
the entanglement entropy of a Fermi liquid for arbitrary 
Fermi surface shape and real space geometry. 

This formula is well defined but complex. There are 
a number of simple checks that can be performed. For 
two or more well separated convex regions, the formula 
above reduces to a sum of terms corresponding to the en- 
tanglement entropy of each region taken separately. For 
free fermions with a nested Fermi surface one immedi- 
ately recovers the one dimensional result as all the fiux 
factors become trivial. At finite temperature, all the log- 
arithms In L are replaced by crossover functions of the 
form In (sinh (ttLT/v)) leading to a simple form at high 
temperature that is extensive in the total real space re- 
gion size. This one dimensional crossover form is well 




Vi V2 Vb 

FIG. 1: Geometry of the entanglement entropy for a non- 
convex real space region. The shaded grey region surrounded 
by the thick black line represents a particular non-convex real 
space region. Dotted lines labeled Vi {i — 1,2,3) represent 
particular one dimensional cuts experienced by a patch of the 
Fermi surface with vertical Fermi velocity. Similarly, dashed 
lines labeled Hi {i — 1,2) represent particular one dimen- 
sional cuts experienced by a patch of the Fermi surface with 
horizontal Fermi velocity. Notice how the vertically moving 
fermion switches from an effectively single interval geometry 
in Vi to a two interval geometry in V2 and back to a sin- 
gle interval in V3. On the other hand. Hi and H2 are both 
effectively single interval geometries. 



known for single intervals in one dimensional conformal 
field theory, but it is nontrivial for the multiple interval 
case. It comes from evaluating higher point correlations 
functions of twist fields on a Euclidean torus in the free 
fermion conformal field theory. As before, the simple 
form depends crucially on the fact that we are dealing 
with free fermions; this result does not hold in a general 
conformal field theory. 

The absence of Landau parameters is reasonable given 
the fact that the entanglement entropy (at least for a sin- 
gle interval) doesn't depend on interaction strength even 
for a Luttinger liquid in one dimension. I cannot com- 
pletely rule out the possibility that there is a prefactor 
/ in the crossover function which is a function only of 
LT: f = f{LT) with the hmits f{x -J> 00) ^ 1 and 
f{x -^ 0) ^ constant, but such a prefactor seems un- 
likely. Given the absence of Landau parameters at high 
temperature and the large U{\)°° symmetry, I believe the 
Fermi liquid entanglement entropy is universal as I have 
described. 

There arc also some subtleties in the evaluation of this 
formula. The most basic subtlety is that the number of 
intervals is subject to change depending on the point on 
the real space boundary chosen. These weak singularities 
are nevertheless completely harmless and do not render 
the integral ill defined. It is best to illustrate the formula 
by a simple example; further results of this nature will 








FIG. 2: Multiple interval geometry and flux factors. A se- 
quence of intervals similar to the situation shown in Fig. 1 V2 
with the local Fermi velocity vertical. The dashed grey line 
encloses the real space boundary segments of interest; the 
segments themselves are the black lines which continue into 
dotted lines outside the dashed grey enclosing line. The dark 
grey shaded regions are interior regions of the real space re- 
gion, while arrows at the surfaces of the these regions indicate 
local normals to the real space boundary. 



be presented elsewhere. Take a d dimensional spherically 
symmetric Fermi liquid, and consider a real space region 
consisting of two spheres Ai and A2 of radius R separated 
by a center-to-center distance L A very interesting quan- 
tity is the mutual information between the two spheres 
defined as I{R,l) = Sai + Sa-2 ^ Sai\jA2- This quan- 
tity is quite powerful as it bounds normalized connected 
correlation functions between operators localized in the 
two spheres. We now evaluate this quantity for the two 
sphere geometry in the limit £ ^ R thus bounding the 



long distance decay of correlations in the Fermi liquid. 

Because we subtract the entropy of each region sepa- 
rately, the only contribution to the mutual information 
comes from configurations that involve both spheres. As 
shown in Fig. 2, each pair consisting of a Fermi surface 
point and a point on the real space boundary only con- 
tributes if the line drawn from the real space point par- 
allel to the local Fermi velocity intersects both spheres. 
The fraction of pairs of points satisfying this constraint 
vanishes in the limit d/ R -^ cxd which is the statement 
that the mutual information vanishes for infinite sphere 
separation. Our task is to estimate how this fraction 
vanishes in said limit. Let one sphere sit at z = and 
the other sit ai z — l\ both spheres sit aX Xi — {) for 
i = l...d — 1 with Xd — z. Since I have assumed spher- 
ical symmetry, the Fermi surface is sphere of radius kp 
in momentum space with kp fixed by density. Now the 
Fermi surface points with vp — ±vpz always gives lines 
connecting the two spheres, but as £/R -^ oo, only a 
small neighborhood around these points continue to gen- 
erate lines connecting the two spheres. We can estimate 
the linear size Sk of this neighborhood as Sk ^ kp {^j^. 
For i/ R — > (X), the difference of logarithms appearing the 
mutual information nearly cancels with a residual (R/l)'^ 
falloff. Putting everything together, I find that roughly 
{RSk^^^ patches on the Fermi surface each contribute 
(R/l)'^ to the mutual information for a total mutual in- 
formation going like 
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As I have already indicated, the mutual information 
provides a bound on connected correlation functions. 
The precise statement is that the square of the connected 
correlation function of any two operators, normalized by 
the operator norms, is bounded by the mutual informa- 
tion. Thus we learn that connected correlation functions 
of local operators in a Fermi liquid must decay at least 
as X 2- as X — > cxd. Calculations are trivially possible 
for free fermions, and one can directly verify that all two 
point functions of bosonic composite operators do indeed 
decay fast enough to satisfy the bound. The bound is sat- 
urated by the two point function of fermions which decays 

_ d + l 

as X T~ times an oscillating function of kpx. This is not 
completely trivial since the fermion creation operator is 
not strictly speaking a local operator i.e. in the Jordan- 
Wigner representation, but this mild non-locality does 
not to spoil the utility of the mutual information. The 
exact coefficient is calculable by performing the geomet- 
rical integrals outlined above. On the other hand, nested 
Fermi surfaces give a mutual information that only de- 
cays as 1/f^ indicating the possible presence of slowly 
decaying correlations. Indeed, this bound is again sat- 
urated by the fermion two point function which decays 
as 1/x for a nested Fermi surface with x parallel to the 
nesting vector. 

I have also numerically computed the mutual informa- 
tion between two rectangular blocks in a two dimensional 




FIG. 3: Geometry of the mutual information for disconnected 
real space regions. The three one dimensional cuts labeled A, 
B, and C represent three different patches on the Fermi sur- 
face for different angles of the Fermi velocity. Cut A comes 
from a mode with vertical Fermi velocity; it always cuts both 
spheres and hence always contributes to the mutual informa- 
tion in this geometry. Cut B also connects both spheres, at 
least for some choices of real space boundary point. However, 
cut C never intersects both spheres simultaneously and hence 
does not contribute to the mutual information. For the sphere 
geometry shown here, there is always a maximum angle from 
vertical such that cuts beyond that angle never intersect both 
spheres. As described in the text, this maximum angle ap- 
proaches zero as the spheres are taken far apart relative to 
their size. 



free fermion system with a nearly spherical Fermi sur- 
face. The geometrical integrals can be done in the limit 
i? > PF with the result 
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where W is the box width, R is the box length, and 
L is the separation between the boxes. This expression 



agrees with the numerical result to within a few percent 
with the deviation due mostly to corrections in W/R and 
R/L; additionally, the qualitative 1/L^ dependence also 
agrees with the numerical results. 



V. NUMBER FLUCTUATIONS 

Finally, I turn to a description of the number fluctu- 
ations. Here at last Fermi liquid parameters will make 
an appearance. Recall that the basic result on the num- 
ber fluctuations is that they scale as L**^^ In L in a Fermi 
liquid for a real space region of linear size L. For free 
fcrmions the precise prefactor is known both from op- 
erator methods and from the conformal field theory ap- 
proach advocated here. Unlike the result for entangle- 
ment entropy, the number fluctuations do depend on Lan- 
dau parameters as interactions are turned on. To see 
this note that the same connected correlation function 
used to determine the number fluctuations is also related 
to the compressibility at finite temperature. This com- 



pressibility does depend on a Landau parameter in addi- 
tion to the effective mass through the density of states. 
Thus we can already say that the number fluctuations 
AN1{L,T) must depend on the Landau parameters in 
the limit T -^ oo (always keeping T ^ Ep.) In fact, we 
can say much more because in Fermi liquid theory the 
interacting density-density correlation function is related 
to the free correlation function by a series of forward 
scatting diagrams that can be summed to give 1/(1 + Fo) 
where Fq is the familiar i = Landau parameter for 
spin-less fermions. The number fluctuations in a d dimen- 
sional Fermi liquid at zero temperature take the asymp- 
totic form 
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with L the linear size of the convex region considered. 

In the presence of interactions at finite temperature I 
find the familiar crossover form 
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where the In L term has been replaced by its finite T gen- 
eralization In (sinh (ttLT/?;)). As before, one can check 
that in the limit of high temperature one recovers the 
usual number fluctuations proportional to region size and 
depending on Fq and the physical Fermi velocity vp. 

The generalization to non-convex regions or multiple 
regions is also possible in line with the generalization 
of the entanglement entropy discussed above. Indeed, 
the structure of logarithms in pUj) also appears for the 
number fluctuations in a one dimensional Fermi liquid. 
The density-density correlation function is 
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and upon integrating x and y over different intervals I 
find exactly the structure of sums of logarithms appear- 
ing in ([To]). This means that the same geometrical in- 
tegral that gave the multiregion entanglement entropy 
in the higher dimensional Fermi liquid also appears in 
the calculation of the number fiuctuations. Indeed, these 
two quantities only differ in the prefactor of the geomet- 
rical integral describing the interplay of real space and 
Fermi surface geometries. To be precise, in the limit 
of large L with a fixed high energy cutoff, I find that 
AN^I/Sl = -^ i^p ■ In fact, given the essentially iden- 
tical forms of the postulated finite temperature scaling 
functions for entropy and number fluctuations, this ratio 
is actually valid even at flnite temperature. 



VI. DISCUSSION 

I have described the calculation of Renyi entropy, mu- 
tual information, and number fluctuations in Fermi liq- 
uids in any dimension for arbitrary codimension one 
Fermi surface and sub-region geometry. These results 
give a very complete description of the low energy struc- 
ture of quantum information in a Fermi liquid. This basic 
Fermi liquid state which underlies so many materials is 
a truly highly entangled quantum phase of matter. Re- 
markably, the entire structure of quantum information 
in these systems is controlled by a beautiful interplay be- 
tween the geometry of the Fermi surface and the struc- 
ture of one dimensional conformal field theory. 

There is one subtlety which deserves a more careful 
discussion, namely the issue of singularities on the Fermi 
surface. The methods described above apply directly 
to smooth Fermi surfaces, but the formalism is strong 
enough to handle aspects of singular situations as well. 
To give a simple example, consider spinless fermions hop- 
ping on a square lattice at half filling. The Fermi surface 
consists of four straight lines running diagonally between 
the midpoints of the Brillouin zone boundary. There are 
singularities on the Fermi surface at (±7r, 0) and (0, ±7r), 
points where the Fermi velocity vanishes leading to well 
known singularities in the density of states. As an ex- 
ample of the subtleties that may arise, consider the mu- 
tual information between two regions separated by a large 
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distance compared to their size along the y-axis. A ht- 
tle thought reveals that the leading contribution to the 
mutual information from the Fermi surface, in the form 
presented above, vanishes in this geometry. In fact, the 
mutual information is not totally zero; indeed, it cannot 
be since the fermion correlation function doesn't vanish 
in the limit of large y. However, the fermion correlation 
function does decay faster (l/y^) than the result for a 
non-singular Fermi surface (1/y^/^). In this sense, the 
vanishing of the leading "Fermi surface" part of the mu- 
tual information is physical, telling us that there are no 
extremely long ranged correlations such as would come 
from a non-singular Fermi surface. When the Fermi sur- 
face geometry does lead to extremely long ranged cor- 
relations, such as along the nesting vector (vr, tt) in the 
half filled example, then the mutual information formal- 
ism above perfectly captures the requisite decay. The 
simple lesson is that there are always subleading terms 
in all quantities which can become visible if the Fermi 
surface contribution happens to vanish. We must simply 
always remember that the formalism developed in this 
paper captures only the leading Fermi surface contribu- 
tion to the physics. 

It should be clear by now that there are other quan- 
tities of interest that can be calculated for Fermi liquids 
using the one dimensional approach. Particularly inter- 
esting would be to study the behavior of entropy and 
other quantities in a non-equilibrium situation, following 
a quench of some type, for example. Powerful techniques 
from one dimensional conformal field should give a great 
deal of control over the non-equilibrium structure of en- 
tanglement in Fermi liquids. There is still the question 
of experimental observation of the effects described here. 
The number fiuctuations are in principle easy to observe. 



we simply need to prepare a clean Fermi liquid, perhaps 
in an optical lattice, and repeatedly count the number 
of fermions in a given subsystem to evaluate AN'^. If 
the subregion can be made large enough I have a quite 
precise prediction for AN'^. Another route is to actu- 
ally physically affect the separation of a given subregion 
from a larger system and measure the resulting entropy 
thermodynamically. I am in the process of studying the 
details of this process. 

Given their experimental ubiquity, it is gratifying to 
have control of the quantum information content of Fermi 
liquids. Of course, in the case of Fermi liquids in the solid 
state, the logarithms will be cutoff by the mean free path 
in the presence of disorder. Still, any clean Fermi liq- 
uid system with a long mean free path will display the 
anomalous quantum information theoretic properties de- 
scribed here. In the larger framework of many-body en- 
tanglement, Fermi liquids have taught us how to produce 
one class of highly entangled quantum states. However, 
much remains to be understood about the role of en- 
tanglement in deeply quantum mechanical phases, espe- 
cially gapless phases like those explored here^Si^"— . Like 
Fermi liquids, these phases should violate the boundary 
law for entanglement entropy. And the search continues 
for other highly entangled phases of quantum matter. 
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